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Abstract 

One-dimensional asymmetric simple exclusion processes (ASEPs) which are coupled to external 
reservoirs via diffusive transport are studied. These ASEPs consist of active compartments char- 
acterized by directed movements of the particles and diffusive compartments in which the particles 
undergo unbiased diffusion. Phase diagrams are obtained by a self-consistent mean field approach 
and by Monte Carlo simulations. The diffusive compartments act as diffusive bottlenecks if the 
velocity of the driven compartments or ASEPs is sufficiently large. A diffusive bottleneck at the 
boundary of the system leads to the absence of a maximal current phase, while a diffusive bottle- 
neck in the interior of the system leads to a new phase characterized by different densities in the 
two active compartments adjacent to the diffusive one and to a maximal current defined by the 
bottleneck. 
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I. INTRODUCTION 



Asymmetric simple exclusion processes (ASEPs) are simple one-dimensional driven lattice 
gases with hard core exclusion. They were originally introduced in the context of protein 
synthesis jl| and have attracted much interest during the last years as simple models for 
wundary-induced phase transitions |2j , for which many rigorous results have been obtained 
•J LlLjJ In the open system, different stationary states are found, which depend on the rates 
of injection and extraction of particles at the ends. Varying the injection and extraction 
rates, or equivalently the densities at the left and right boundary, both continuous and 
discontinuous phase transitions are observed. The actual stationary state is selected via the 
dynamics of domain walls and density fluctuations J(|. 

Promising candidates for the experimental observation of these phase transitions are 
systems of cytoskeletal motors which move unidirectionally along cytoskeletal filaments [?], 
0, |^ . However, these motors unbind from their track after a few seconds, since their binding 
energy can be overcome by thermal fluctuations. Observed over sufficiently long times which 
exceed a few seconds, they alternate between the bound and the unbound state and perform 
peculiar random walks. If a motor is bound to a filament, it moves in a directed way along 
the filament, while unbound motors diffuse freely. As motors are strongly attracted by 
the filament, the motor density along the filament can be large even if the overall motor 
concentration is rather small, which implies that hard core exclusion plays an important 
role in the bound state. To study these combined movements, we have recently introduced a 
class of lattice models where bound and unbound motor movements are described as biased 
and symmetric random walks on a lattice, respectively 

BBS 

. In these models, the traffic 

of motors along a filament is an asymmetric simple exclusion process with the additional 
property that motors can attach to and detach from the track. 

For open tube systems with a single filament and fixed motor concentrations at the tube 
ends, the same types of phases are found as for the usual one-dimensional ASEP j^. If the 
filament is shorter than the tube and motors have to diffuse over a certain distance to reach 
one end of the filament from the left boundary and again to reach the right boundary from 
the other end of the filament, the phase boundaries of the system can be shifted by changing 
geometrical tube parameters or motor parameters. In particular, a maximal current phase, 
in which the current attains its maximally possible value, can only be found if the diffusive 
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currents from the left end of the tube to the filament and from the filament to the right end 
of the tube can be as large as this maximal current. These diffusive currents however are 
restricted by the diffusion coefficient of the unbound motors and by geometric parameters 

D 

The latter phenomenon is generic and not restricted to the specific tube geometry. In the 
following, we study several one-dimensional systems which consist of compartments char- 
acterized by active or diffusive transport. We will discuss four simple geometries as shown 
schematically in Fig. ^ While particles move only to the right in the active compartments, 
forward and backward steps occur with the same probability in the diffusive compartments. 
For these models, we determine the phase diagram analytically using a mean field approach 
to calculate effective boundary densities or effective injection and extraction rates for the 
active compartments. The method is based on the constraint that the stationary current 
must be equal in all compartments. 

The article is organized as follows. After introducing the model in section |nj we discuss 
diffusive injection and extraction of particles into/from an active compartment in section 
IIHI We start with only diffusive injection in section IlII Al which corresponds to case (A) in 
Fig. m proceed with only diffusive extraction in section IIIIBI see case (B) in Fig. Q and 
then study the case (C), for which particles are both injected and extracted via diffusive 
compartments, see section IIII CI We compare the mean field results with Monte Carlo 
simulation in section IIII Dl Finally, we discuss the case of a diffusive compartment between 
two active compartments as shown as case (D) in Fig. ^ in section HVl This case corresponds 
to a defect that must be overcome by unbiased diffusion. 

II. THE MODEL 

In the following, we will discuss transport on one-dimensional lattices. The coordinate 
along the lattice is denoted by x and will be measured in units of the lattice constant i. 

We consider systems that can be decomposed into two or three different compartments 
in which transport is either diffusive or directed. The four cases that will be discussed in 
the following are shown schematically in Fig. ^ The total length of the system is taken to 
be L in all cases. The linear extensions of the compartments are denoted by L 1; L 2 , and L 3 , 
compare Fig. ^ In case (A), the system consists of a left compartment with 1 < x < L\, 
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where transport is diffusive, and a right compartment with Li + 1 < x < L\ + L 2 = L, where 
transport is active or directed. In case (B), transport is directed in the left compartment 
and diffusive in the right compartment. The cases (C) and (D) correspond to situations 
where the systems consist of three compartments with L = L\ + L 2 + L 3 . In case (C), 
transport is driven in the middle compartment with L\ + 1 < x < L\ + L 2 and diffusive in 
the left and right ones. Finally in case (D), transport is directed in two compartments, the 
left and the right one, but diffusive in the middle compartment. In all cases, we will assume 
that the extensions of the active compartments are sufficiently large, so that we can neglect 
finite-size effects. 

In the following, we will take the active transport to be always directed to the right and 
to be totally asymmetric, i.e., we do not allow backward steps in the compartments with 
active transport. At lattice sites which belong to such an active compartment, particles 
attempt to hop to the adjacent lattice site to their right with a certain probability per unit 
time r. We denote this probability by v since it is equal to the velocity of a particle in the 
active compartment (and in the absence of other particles), measured in units of £/t. The 
hopping attempt is rejected if the target site is already occupied by another particle. In 
summary, motion in the active compartment is described by a totally asymmetric simple 
exclusion process. 

In contrast, motion in the diffusive compartments is described by a symmetric exclusion 
process. A particle at a lattice site which belongs to a diffusive compartment attempts to 
make a forward and a backward step with equal probability D, which corresponds to the 
diffusion coefficient measured in units of £ 2 /t. Note that we could eliminate one parameter 
by measuring time in units of the time scale for diffusive steps of size £ by choosing r = 
£/(2D) (for this choice, the diffusion coefficient measured in units of £ 2 /r would be given 
by D = 1/2). This implies that the results which we derive in the following will depend 
only on the ratio v/D. All hopping attempts in the diffusive compartments are again 
rejected if the target site is occupied by another particle. In order to simplify the following 
calculations, we do not allow particles to enter the active compartments from the right, i.e. 
all hopping attempts from the first site of a diffusive compartment to the last site of the 
active compartment to its right — from L x + 1 to L x in case (B) and (D) and from L x + L 2 + 1 
to Li + L 2 in case (C) — are rejected. 
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Finally, the densities at the boundary sites x = and x = L + 1 have the fixed values 

p(x = 0) = and p(x = L + 1) = p ex . (1) 

These sites are taken to have the same dynamics as the adjacent compartments of the system. 
Particles thus attempt to enter the system from the left with probability Dp m if the first 
compartment is diffusive and with probability vp- in if it is an active compartment. Particles 
at the last lattice site with x = L leave the system to the right with probability v(l — p ex ) 
if the site x = L belongs to an active compartment and with probability D(l — p ox ) if it 
belongs to a diffusive one. In the latter case, particles also try to enter the system from the 
right at x = L with probability Dp cx . Likewise, particles at site x — 1 can leave the system 
with probability D(l — p m ) if x = 1 belongs to a diffusive compartment. As before, particles 
can only enter the system at x = 1 or x = L if these sites are not occupied. 



III. DIFFUSIVE INJECTION AND EXTRACTION OF PARTICLES 

In this section, we consider the cases (A)-(C), where transport is driven or active in one 
compartment, but particles are diffusively injected and/or extracted from the system and 
have to diffuse over a certain distance before they reach the active compartment and/or 
before they can leave the system at the right end. 

As mentioned before, the active compartment is described by an asymmetric simple 
exclusion process. Let us therefore briefly summarize what is known about the phase diagram 
of this process, see, e.g., Ref. Q. In an open system, where the densities are fixed to p m 
and p cx at the left and right boundary of the ASEP, respectively, the stationary state is 
determined by the boundary densities. The stationary states are characterized by the bulk 
density p° and the stationary current J. 

For the ASEP with open boundaries, three different phases can be distinguished as shown 
in Fig. EJ If the density at the left boundary is relatively small and satisfies p in < 1/2 and 
Pi n < 1 — p ex , the system is in the low density phase (LD) for which the bulk density is 
equal to the left boundary density and the current is J = t>pi n (l — Pin)- If the density at 
the right boundary is relatively large with p cx > 1/2 and p cx > 1 — p in , the system is in 
the high density phase (HD), the bulk density is equal to the right boundary density, and 
the current is J = vp cx (l — p cx ). At the transition from the low density to the high density 
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phase, the current is continuous, but the bulk density is discontinuous. Finally, for p- m > 1/2 
and p ex < 1/2, the system is in the maximal current phase, where the current is maximal, 
J = v/4 and the bulk density is 1/2. The transitions to the maximal current phase are 
continuous. 

In contrast to the simple ASEPs just described, our systems are characterized by the 
property that at least one of the boundary densities of the active compartments is not fixed, 
but adjusted by the dynamics of the system. In the following, we will determine the phase 
diagrams of these systems using a mean field approach. We proceed in three steps and 
consider first diffusive injection and extraction of particles separately, combining them in 
the last step. 



A. Diffusive injection of particles 

First we consider case (A), a system with only diffusive injection of particles. Particles 
leave the system at the right boundary with rate v (1 — p ex ) and no particles enter the 
system at the right end. In the stationary state, the current J must be the same in both 
compartments. In the left compartment with 1 < x < L% where transport is diffusive, the 
density is then given by p(x) = p m — xJ/D. Within mean field approximation, the right 
compartment with L\ + 1 < x < L, corresponds to the usual ASEP with the effective left 
boundary density 

_ Dp(Lx) _ Dp in LiJ 

Pin,eff — — • {*) 

V V V 

as follows from vp- infi ^ = Dp(L{). The quantity upi n , e ff corresponds to the rate with which 
particles attempt to enter the ASEP at its left boundary. Note that (i) this effective bound- 
ary density depends on the current J and (ii) that pi n , e ff can be larger than one. The right 
boundary density is given by p cx . 

The phase diagram can now be determined in a self-consistent way. As in the tube 
system studied in Ref. the same phases are found as for the one-dimensional ASEP, but 
the location of the transition lines depends on the values of the model parameters v/D and 
Li, and the maximal current phase may be shifted out of the physically accessible range of 
the parameters. 

The system is in the maximal current phase if p ex < 1/2 and pi n ,eff > 1/2. In this case 
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the current is J = v/4, and the condition 



D Lx 1 

Pin,eff = —Pin - -J- > ^ I 3 ) 



implies that the system is in the maximal current phase for 
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For large f/-D, the latter value of the left boudary density is larger than one and therefore 
not physically accessible. This implies that a maximal current phase is only present for small 
velocities with v/D < 2/(l+Li/2). If the velocity is larger, unphysically high densities would 
be necessary at the left boundary in order to establish a sufficiently large density gradient 
which could generate a diffusive current with the value v/4, the maximal current defined by 
the driven compartment. In this situation, the diffusive compartments acts as a diffusive 
bottleneck: If the maximally possible diffusive current through the diffusive compartment is 
smaller than v/4, a maximal current phase cannot occur, because the diffusive compartment 
cannot maintain the maximal current. 

A simpler estimate comparing the maximal diffusive current D/L\, which is restricted 
by the maximal density difference of one, with the maximal driven current v/4 yields the 
condition, v/D < 4/Li, which agrees with the previous one for large L\, but is less restrictive 
for small L\. The latter discrepancy reflects the fact that the maximal density difference in 
the diffusive compartment is actually smaller than one since the density at x — L\ must be 
larger than zero. 

In addition, a low density phase is found for Pi niC ff < 1/2 and Pi n , e ff < 1 — p cx and a high 
density phase for pi n>c fr > 1 — p cx and p ex > 1/2. Along the transition line between the high 
density and low density phases we can use J = fp ex (l — p ex ) and obtain 



1 + (L, - l)p cx - L lP l (5) 



for the transition line between the low density and the high density phase. This line extends 
from the right upper corner of the phase diagram to the right upper corner of the maximal 
current phase region. If there is no maximal current phase, it ends at a point with p in = 1 
and p ex > 1/2. 

Phase diagrams for two cases are shown in Fig. El We have chosen L\ = 10 and D = 1/2 
in both cases. The condition for the presence of a maximal current phase is then v < 1/6. In 



Fig. Efa), the velocity is v = 0.1 < 1/6 and all three phases are present, while in Fig. Efb), 
v = 0.2 > 1/6 and the maximal current phase is absent. In the latter case the largest part 
of the phase diagram is covered by the low density phase. 

In the maximal current phase the current is J — v/A. In the high density phase, it is 
determined by the right boundary density p ex and has the value J = t>p ox (l — p ox ). Finally 
in the low density phase, the current is given by the self-consistency condition 

J = Vp in , eS (j)[l - p inj eff (J)], (6) 

which leads to a quadratic equation for the current. The solutions is uniquely determined 
by the limits J = for p m = and J = t>/4 for p in = ^(l + L\/2) and is given by 



B. Diffusive extraction of particles 

Next we consider case (B), in which particles which reach the end of the active compart- 
ment have to diffuse over a distance L 2 before they can leave the system at the right end. 
This case can be treated in the same way as the one with diffusive injection. Note, however, 
that it cannot simply be derived from the the latter one using particle-hole symmetry, be- 
cause a particle at the site left of the driven compartment attempts to enter it with rate D, 
while a hole at the site right of the driven compartment does so with rate v. 

The density profile in the diffusive compartment is given by p(x) = p cx + (L + 1 — x)J/D 
for Li + 1 < x < Li + L 2 = L. Therefore the effective rate with which particles attempt to 
leave the active compartment at x = L\ is v (1 — p cx — L 2 J / D) = v(l — p cx , c g) corresponding 
to an effective right boundary density of the driven compartment given by 

Pex,efT = Pex + L 2 J/D. (8) 

The maximal current phase is now found for p in > 1 /2 and 

1 L 2 v , N 

Pex < - - W , (9) 

which is always < 1/2. Again the maximal current phase is only present if the range 
of boundary densities defined by Eq. Q is physically accessible. Here the corresponding 
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condition is p cx > 0, which is valid for v/D < 2/L 2 . The latter condition expresses again the 
fact that the diffusive compartment must also support this maximal current. The diffusive 
current is, however, restricted by the maximally possible value of the density gradient in 
the right (diffusive) compartment, 1/(2L 2 ), which leads to a maximal diffusive current of 
D/(2L 2 ). If the latter current is smaller than v/4, a stationary maximal current phase is 
absent; therefore, the presence of the maximal current phase requires that the maximal 
diffusive current is > v/A which leads to the condition v/D < 2/L 2 . 

The condition 1 — p cx , e fr = Pin with J = t>Pi n (l — Pi n ) yields the transition line between 
the low density and the high density phase 

Pcx = 1 " Pin(l + ^L 2 ) + ^L 2 pf Q . (10) 

For velocities larger than 2D/L 2 , which is the maximal value for the occurrence of a maximal 
current phase, this line ends at a point in the phase diagram with p in = and p ex < 1/2. In 
this case the high density phase covers most of the phase diagram. 

The current is J = v/A in the maximal current phase and J = vp in (l — p in ) in the low 
density phase. In the high density phase, it is again given by a self-consistency condition 
J = "ypex,eflf(l — Pex.cfr), where pex.efr is a function of J, from which we obtain 



J 



(V / V V /f\ 2 \ , 

-l + -L 2 (l-2p cx ) + yi-2-L 2 + 4-L 2 p cx +(^-L 2 j 1. (11) 



2(L 2 v/Df 

C. Both diffusive injection and extraction of particles 

Now we consider case (C), i.e., we combine the two preceding cases. Transport is now 
driven in the middle compartment for which we have the effective boundary densities 

Pin,eff = - — and p CXjCfr = p cx + L 3 J/D. (12) 

v v 

The phase boundary between the low density phase and the maximal current phase is not 
affected by adding another compartment at the right end, thus we can use the result from 
case (A). Likewise the phase boundary between the high density phase and the maximal 
current phase is unaffected by the left diffusive compartment, so that we can use the result 
from case (B) upon substituting L 2 with L 3 . The maximal current phase is therefore found 
for 

v ( L\\ 1 L 3 v 

P^>2d{ 1 + Y) and Pcx< 2"4D- (13) 
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It is present if v/D < 2/(1 + L\/2) and v/D < 2/L3. In addition, the current is given by 
Eqs. (J7J) and (fTTj) in the low density and the high density phase, respectively. 

Finally, the transition line between these two phases is obtained from the condition 

Pin,eff = 1 - Pcx,cfr, which leads to 



where J(p 

ex) is the current in the high density phase for the right boundary density p ex as 
given by Eq. (fTT]) . The complete phase diagram is shown in Fig. where we have chosen 
parameters for which a maximal current phase is present. 

D. Comparison with simulations 

In addition to the self-consistency calculations presented above, we performed Monte 
Carlo simulations. In this section, we compare the simulation results with the mean field 
predictions for case (C). 

In the case where the maximal current phase is absent, i.e. for large velocities, we find 
quantitative agreement of the measured current and bulk density with the predictions of 
the mean field calculation. The transition from the low density to the high density phase 
occurs at the predicted density. As an example, we show results for the bulk density p\ in 
the active compartments in Fig. where we have chosen p ex = and p ex = 0.5. In the 
first case, the system is in the low density phase for all values of p in , in the second case, a 
transition to the high density phase is found at p m ~ 0.54. The mean field results (lines) 
and the simulation data (symbols) agree well. 

If a maximal current phase is present, i.e. for small velocities, the agreement is less good, 
although the phase diagram is still in qualitative agreement with the mean field predictions. 
Fig. El shows again results for the case p ex = 0. Close to the transition to the maximal 
current phase the current is smaller than predicted by the mean field calculation. Therefore 
the transition to the maximal current phase is shifted towards a larger value of p m and 
the increase of the bulk density near the transition is less steep. Far from the transition, 
however, agreement is again good. Likewise, the transition line between the low density and 
the high density phase is also shifted towards larg GT Pin; this transition line ends on the 
phase boundary of the maximal current phase. Agreement becomes again better far from 




(14) 
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the maximal current phase, since the other end point of the line (p- m = 0, p ex = 1) is exact. 



IV. DIFFUSIVE BOTTLENECK IN THE MIDDLE 

Finally, let us consider case (D), where active transport is interrupted by a diffusive 
compartment. In the case of molecular motors, this can be realized by a gap in the filament 
network, along which active transport takes place. Motors thus have to overcome this gap 
by diffusion before they can rebind to a filament and continue their active movements. If 
the middle compartment consists of only one lattice site, L 2 = 1, this system reduces to the 
case of an ASEP with a single defect which has been discussed previously, see Refs. ldL lll|. 

The effective right boundary density for the left active compartment is 

Pex,efT = p(Ll + 1) (15) 

and the effective left boundary density for the right active compartment is 

DpjL, + L 2 ) Dp(L 1 + l)-(L 2 -l)J 

Pin,cfT — — ■ l-LDJ 

V V 

In this case, there are five possible phases. Because the current is the same in both active 
compartments, the bulk densities in the left and right compartment, p° and p°, respectively, 
are either equal or related by p\ = 1 — If the bulk densities in both active compartments 
are equal, p\ = p[j, there are three possibilities. Both compartments can be in the low 
density, high density or maximal current phase. We denote these three cases by LD-LD, 
HD-HD, and MC-MC, respectively. If the densities are not equal, we have p\ — 1 — p°, and 
there are two additional possible phases, where one compartment is in the high density and 
the other in the low density phase. These phases will be denoted by HD-LD and LD-HD if 
p\ is larger or smaller than p®, respectively. 



A. Phases with equal densities in the active compartments 

In the MC-MC phase, the current is J = v/A and we have four conditions for the boundary 
densities, p in > 1/2, p cx < 1/2, pi n , e fr > 1/2 and p ex , e s < 1/2. The first two conditions are 
the same as for an ASEP without the diffusive compartment and the latter two conditions 
yield v/D < 2/(1 + L 2 ). The maximal current phase is found only for small velocities, for 
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larger velocities in the active compartments, the diffusive section acts again as a diffusive 
bottleneck. 

The LD-LD phase is characterized by J = vp in (l — p in ) = t>Pm, e ff(l ~ Pm, e ff)- As both p in 
and Pin je fT must be smaller than 1/2, this implies p m = pi n , e ff- Together with the condition 
p ex < 1 — Pin, e ff for the right active compartment, this implies that the LD-LD phase is found 
within the region of the phase diagram where the low density phase of the usual ASEP is 
located. An additional condition is given by p ex ,ofr = p{L\ + 1) < 1 — p in . Thus, with the 
condition p in = Pi n;0 ff, we obtain p{L\ + 1) as a function of p in , 

p(L x + 1) = ^[p in + (L 2 - l)p in (l - p in )] (17) 

which leads to the inequality 

"^[Pin + (L 2 - l)p in (l - Pin)] + Pin < 1. (18) 

Because p m < 1/2, the left hand side of this inequality is increasing monotonically and the 
solution is given by p; n < pi n ,*, where p- in ^ is the solution of the equation which is obtained 
upon substituting '<' with '=' in Eq. (|18jl. This solution is uniquely determined by the 
limiting case L 2 = 1, in which Eq. (J18)) yields p; n < 1/(1 + v/D). As a result we obtain the 
condition 



^ 2 + l- y/(^ 2 + l) 2 - 4%(L 2 -1) 

Pin < Pin, - 2 ^(L 2 - 1) • (W) 

If pi nj * is larger than 1/2, this condition does not further restrict the occurrence of the 
LD-LD phase, since for p in = 1/2 the transition to the maximal current phase takes place. 
Indeed, pi n) * > 1/2 implies the simpler condition v/D < 2/(L 2 + 1), which is exactly the 
condition we derived above as a condition for the presence of the MC-MC phase. 

On the other hand, if p- in ^ < 1/2, condition ()18|) yields a restriction of the LD-LD phase 
to the region in the phase diagram with p; n < p- m ^. Therefore we conclude that for large 
velocities with v/D > 2/(L 2 + 1), the transition to the MC-MC phase at p; n = 1/2 is 
replaced by a transition to another phase at p; n = pin,*- The only possibility for this phase 
is the HD-LD phase which will be discussed below. 

Finally note that for L 2 = 1 our result recovers the condition for the ASEP with a defect 



11 1, where for v/D < 1 the phase diagram is the same as without the defect, while for 
v/D > 1 a phase dominated by the defect is found. 
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The HD-HD phase can be treated in the same way and gives corresponding results. A 
HD-HD phase can be found for p ox > 1 — p m and p ex > 1/2 for small velocities which fulfill 
again the condition v/D < 2/(L 2 + 1), while for larger velocities an additional restricting 
condition is found, namely p cx > p eXj * with p eXJk = 1 — p m ^. 

B. Phases with different densities in the active compartments 

For an LD-HD phase, the current must be J = t>pi n (l — pin) = t>p ex (l — Pcx)- This implies 
that the densities are p in = 1 — p ex , i.e., a stationary state with a low density in the left, but a 
high density in the right active compartment can only be expected along the line, where the 
low density and high density phases coexist in the ASEP without a diffusive compartment. 
In this case, however, a domain wall diffuses through the system in the usual ASEP resulting 
in a density profile which increases linearly [6j]. The same can be expected for our case with 
the exception of the regions close to the diffusive compartment. This behavior has previously 
been observed in simulations for the case of an ASEP with a defect ^j] which correspond to 
our case with L 2 = 1 and we find the same behavior in Monte Carlo simulations for larger 
L 2 , see Fig. 

Finally, let us consider the possibility of a HD-LD phase. In this case the current is 
J = vp in>eS (l - Pin.cfr) = vp CX)CS (l - p cx ,eff) which, together with the conditions p injCfr < 1/2 
and p ex , e ff > 1/2, implies pi n , e ff = 1 — Pex, e ff- Substituting the expression for the effective 
boundary densities, we obtain a quadratic equation for p{L\ + 1) with the solution 

p(Li + 1) = Pcx,* = 1 - Pin,* (20) 

with pi ni!t as given by Eq. (|19p. For L 2 = 1 we recover again the corresponding result for the 
ASEP with a defect. 

C. Phase diagrams 

We can now summarize the results into phase diagrams as shown in Fig. |HJ There are two 
different cases: If v/D < 2/(1 + L 2 ), the phase diagram corresponds to the one of the ASEP 
without the diffusive compartment, see Fig. [S^a). The bulk density is equal in both active 
compartments. If, on the other hand, the velocity is larger with v/D > 2/(1 + L 2 ), the 
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diffusive compartment acts as a bottleneck and the phase diagram is modified, see Fig. |Sfb). 
We find the LD-LD phase for p in < p in ^ and p in < 1 — p ex and the HD-HD phase for 
Pex > Pex,* and p cx > 1 — p in . For p in > p in * and p ex < p CXj *, the system exhibits the 
HD-LD phase. While in the LD-LD and HD-HD phases, the bulk densities in the active 
compartments and the current are determined by the boundary densities, these quantities 
are independent of the boundary densities in the HD-LD phase and depend only on the 
ratio v/D and the length L2 of the diffusive compartment. 

The HD-LD phase has some similarities to a maximal current phase. The current is 
constant throughout this phase and attains its maximal value compared to all other phases, 
J = v Pm,*(l — Pin,*) = ' u Pcx,*(l _ Pcx,*)- This maximal value, however, is not determined by 
the active compartments, but corresponds to the maximal current which can be supported 
by the diffusive compartment. Correspondingly, the density profiles in the HD-LD phase 
(shown in Fig. EJ) do not exhibit the power-law behavior known from the usual maximal 
current phase. The bulk densities in the left and right active compartments are constant in 
this phase as well and are given by p ox ,* and pi n ,*, respectively. Note, however, the following 
difference compared to the usual maximal current phase. The transitions to the maximal 
current phase in the usual ASEP are continuous. The transitions to the HD-LD phase in 
our case are somewhat peculiar in the sense that they are continuous in one compartment, 
but discontinuous in the other. For example, at the transition from the LD-LD phase to the 
HD-LD phase, the density p° in the right active compartment is continuous, but the bulk 
density p° in the left active compartment exhibits a jump from p in ^ to 1 — pi n ,*- 

We performed again Monte Carlo simulations and compared the results with the mean 
field predictions. As a result, we find that the phase diagram obtained by the self-consistent 
mean field approach is recovered for small velocities. For large velocities, on the other 
hand, the qualitative behavior is correctly predicted by the mean field calculation, but the 
transition lines for the transition to the HD-LD phase are shifted. For example for v = 1, 
D = 1/2, and L2 = 3, we find the transition from the LD-LD phase to the HD-LD phase 
at pi n ~ 0.08 in the simulations, while our mean field calculation predicts a transition at 
Pm = Pin,* — 0.13. Correspondingly, there are also differences in the values for the current in 
the HD-LD phase and the critical value of the velocity is found to be smaller than predicted 
by the mean field calculation. For L 2 = 3 and D = 1/2, we observe the usual maximal 
current phase for v ^ 0.15 in the simulations while the mean field approach yields v < 0.25. 
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Finally, let us add a remark concerning the density profiles in the HD-LD phase. From 
our mean field approach, we expect the constant bulk densities in the left and right active 
compartments to be approached exponentially from the left and right boundary, respectively 
This is the case in the simulation data; however, in addition, an excess density close to the 
diffusive compartment is observed, which is not expected from the mean field approach, see 
Fig. El As reported previously for the case of an ASEP with a single defect site , this 
excess density decays as a power-law ~ x~ x , so that some kind of long-range order is also 
present in this phase which plays the role of a maximal current phase for the system with a 
diffusive bottleneck. 



V. SUMMARY AND CONCLUSIONS 



We have discussed transport in one-dimensional lattices which consist of two or three 
compartments where transport is either driven or diffusive, a situation that is inspired by 
the motion of molecular motors which diffuse until they reach a filament and then move along 
that filament in a directed manner Mutual exclusion from lattice sites is important, as 
many particles can be injected into the system from reservoirs of fixed density at both ends. 
This is again realistic for molecular motors, which are strongly attracted by the filaments, 
so that the density of motors along the filaments can become large, even if the motor 
concentration in solution is relatively small. Traffic in the compartments where transport is 
active or driven is described by an asymmetric simple exclusion process. 

We have studied four different geometries. In the cases (A)-(C), particles are injected 
into an active compartment and/or extracted from it via diffusive compartments. In case 
(D), active transport is interrupted by a diffusive compartment in the middle. The latter 
case is a generalization of the ASEP with a point defect. 

In all cases, the diffusive compartments can act as diffusive bottlenecks. If the velocity of 
the particles in the driven compartment is sufficiently small, the phase diagram is essentially 
the same as for the usual one- dimensional ASEP. In the cases (A)-(C), the locations of the 
transition lines depend on the model parameters. For large velocities, on the other hand, 
transport through the lattice is limited by the diffusive compartments which cannot support 
arbitrarily large currents. In the cases (A)-(C), this situation implies the absence of the 
maximal current phase. 
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In contrast, in case (D), this leads to a peculiar new phase, the HD-LD phase, where the 
density is high in the left, but low in the right active compartment. As in a maximal current 
phase, the current is constant in the HD-LD phase and has the maximal possible value. In 
contrast to the usual maximal current phase, this value of the currrent is determined by the 
diffusive compartment, and the transitions to the HD-LD phase are continuous in one, but 
discontinuous in the other active compartment. 
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FIG. 1: Four geometries of one-dimensional lattices which consist of active compartments (gray) 
where transport is described by an asymmetric simple exclusion process and diffusive compartments 
(white) where transport is described by a symmetric exclusion process. The linear extensions of 
the compartments are denoted by Li, L2, and L3. The total length of the system is given by 
L = L\ + L2 in cases (A) and (B) and by L = L\ + L2 + L% in cases (C) and (D). In the active 
compartments, motion is completely biased and particles hop to the right with probability v, while 
in the diffusive compartments, particles hop both to the left and to the right with probability D. In 
addition, we do not allow particles to enter the active compartments from the right. The densities 
at the boundary sites x = and x = L + 1 are fixed to p- m and p cx , respectively. 



17 





Pin 

FIG. 2: Phase diagram of the usual asymmetric simple exclusion process (ASEP), which describes 
a single active compartment, as a function of the left and right boundary densities p m and p ex . 
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FIG. 3: Phase diagrams for the ASEP with diffusive injection of particles at the left boundary, i.e. 
for case (A) as shown in Fig. ^ as a function of the left and right boundary densities pi n and p cx . 
The parameters are (a) v = 0.1, D = 1/2, and L\ = 10; and (b) v = 0.2, D = 1/2, and L\ = 10. 
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FIG. 4: Phase diagram for the ASEP with diffusive injection and extraction of particles, i.e. for 
the case (C) shown in Fig. ^ The parameters are v = 0.1, D = 1/2, and L\ = L-j, = 6. 
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FIG. 5: Bulk density p\ in the active compartment of a system with geometry (C) as a function 
of the left boundary density p in . Lines are the mean field results and symbols are simulation data. 
The right boundary density is fixed to p ex = (solid line) and /0 ex = 0.5 (dashed line). Parameters 
have been chosen so that no maximal current phase is found; v = 1, D = 1/2, L\ = L3 = 6 and 
L 2 = 388. 
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FIG. 6: Current J (dashed line) and bulk density p\ (solid line) in the active compartment of a 
system with geometry (C) as a function of the left boundary density p- m for a case with transition 
to the maximal current phase. The lines are the predictions of the mean field calculation, the 
symbols simulation data. The parameters are v = 0.1, D = 1/2, L\ = L3 = 6, L2 = 388, and 
Pcx = 0. 
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FIG. 7: Density profile for the 'LD-HD' phase or coexistence of the LD-LD and HD-HD phases 
for a system with a diffusive defect [case (D)]. As discussed in the text, the profile is essentially 
linear rather than consisting of two plateaus. The parameters are v = 1, D = 1/2, L = 205, 
L\ = L3 = 101, L2 = 3, pi n = 0.05, and p cx = 0.95. 
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FIG. 8: Phase diagrams for case (D), i.e. for a diffusive compartment between two active compart- 
ments, as a function of the left and right boundary density, p m and p ex . (a) Small velocity, v = 0.1 
and (b) large velocity, v = 1. The diffusive section has length L2 = 3 and the diffusion coefficient 
is D = 1/2. 
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FIG. 9: Density profile for the HD-LD phase in case (D) with p in = 0.5, p cx = 0, v = 0.25, and 
D = 1/2. The geometrical parameters are as in Fig. [7| 
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